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I. INTRODUCTION 

The ATLAS [1] and CMS [2] experiments at the Large Hadron Colhder (LHC) have recently 
presented results that indicate the observation (also supported by evidence coming from the 1.96 
TeV run of the Tevatron [3j) of a new resonance. When interpreted as a Standard Model (SM) 
Higgs, the combined channels are consistent with the SM expectation at la, yet individual channels 
show deviations from the SM expectation in the l-2a range. The discovery of the Higgs boson 
would certainly be one of the most exciting developments in particle physics to date and it is 
tempting assume this new resonance is indeed the Higgs, but establishing the true nature of this 
excess as the SM Higgs must still proceed with due diligence. 

The two main theoretical inputs in performing such a Higgs identification are the Standard 
Model branching fractions for each of the decay channels used in the combination and the overall 
Higgs production cross section. We highlight that Higgs production from gluon fusion, the dom- 
inant production mode at hadron colliders, occurs via loops of SM quarks and is hence uniquely 
sensitive to New Physics (NP) effects arising from new colored states [1H26] or more general Higgs 
portal |27H39| and scalar mixing effects |40H44| . The 77 branching fraction, which arises at loop- 
level in the Standard Model [45j, is similarly sensitive to NP effects |46H50j . For example, in the 
well-studied four generation Standard Model (SM4), gluon fusion rates are enhanced while the 
branching ratio to diphotons is suppressed 061 [5Tti6T] . 

Moreover, the landmark discovery of a Standard Model Higgs crucially relies on affirming the 
hypothesis of the Higgs mechanism for spontaneous breaking of SU{2)l x U{1)y gauge symmetry, 
which predicts the existence of the Higgs boson. This fundamental test can only come after directly 
measuring the Higgs couplings to the SU{2) gauge bosons. The fact that the Higgs is responsible 
for chiral symmetry breaking in the Standard Model and hence gives fermion masses is only a 
byproduct of the Higgs mechanism, and thus, in particular, the gluon fusion Higgs production 
mode does not directly probe the Higgs mechanism. This implies that NP could still be hiding in 
the gluon fusion process without effecting EWSB. 

This crucial point, which has also been emphasized in several recent papers [I9 l \20 \ [32 | l35l [38] . 
means that the LHC Higgs searches can be skewed by the presence of new colored particles which 
positively or negatively contribute to gluon fusion. Moreover, the excess in the data should be 
interpreted not only in the context of a SM Higgs, but also in the more exciting scenario of a 
possible new scalar state which arises from a Higgs portal-induced mixing between the SM Higgs 
and a new scalar. We demonstrate that extended color sectors involving new colored particles will 
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generally give rise to both effects. In particular, if the new colored particles do not get their mass 
from the vacuum expectation value (vev) of the Higgs boson, then a generic Higgs portal term can 
give rise to Higgs mixing. We can see that direct Higgs coupling and Higgs portal-induced scalar 
mixing are two important categories of NP contributions that can have marked effects on Higgs 
collider signals, and thus we consider them simultaneously. 

Motivated by the possibility of probing new colored states via gluon fusion, we adopt a building 
block approach for an arbitrary NP model. Namely, we isolate and calculate the gluon fusion 
amplitude for new colored scalars, fermions, and vectors. In the case of the colored vector, we 
present the calculation in the context of the renormalizable coloron model (ReCoM) |62H64j such 
that the concomitant effects from maintaining UV consistency can be readily included. We also 
allow for Higgs mixing, where the SM Higgs is mixed with a new scalar. In addition, for a mild and 
well-motivated set of assumptions, we give generic expressions for branching ratios of the scalar 
mass eigenstates into the most sensitive SM Higgs decay modes. 

Since our work has some overlap with many studies in the literature, we survey several repre- 
sentative papers and elaborate on the differences. A few recent papers have focused on our first 
category of NP effects for gg ^ h in which the NP states couple directly to the SM Higgs. In par- 
ticular, the authors of |32j focused solely on the situation where new particle masses arise from the 
Higgs vev, which simultaneously sharpens their discussion of resulting gluon fusion and diphoton 
decay phenomenology and limits the breadth of their conclusions. Separately, the authors of |20j 
focused on Higgs portal phenomenology with new colored scalars, while the work in |33| [35] also 
included new scalars transforming under the full SM gauge symmetry. A similar study, emphasiz- 
ing the constraints from electroweak precision fits, was performed in jl8[ I19j. We go beyond these 
direct coupling studies by also including the effects of a colored vector. 

There have also been a number of recent studies \65\ |66] where fits are performed in order to 
determine how consistent the data is with a SM Higgs hypothesis. These studies find that generally 
the excess is largely consistent with a SM Higgs with a tantalizing, but small enhancement in the 
77 channel. However, these fits still have large uncertainties, and the next data set can change 
the picture drastically. Also, the Higgs couplings to bosons rely heavily on the vector boson 
fusion channel, which has large fiuctuations between 7 and 8 TeV. These uncertainties leave room 
for modifications to the coupling of the SM Higgs to gluons which can be either enhanced or 
suppressed given the sign of the Higgs portal term. We will examine this in detail below. 

Regarding our second category of NP effects, when the SM Higgs mixes with a new scalar via 
a Higgs portal term, a majority of the literature has focused on the case where the New Physics 
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sector is completely invisible to the SM [27H3QI [Ml [Ml EH] , providing a possible connection to the 
dark matter. In this situation, as we will see in Sec. only a simple mixing angle is needed 
to parametrize the effects on Higgs phenomenology, if no new decays are kinematically allowed. 
Our work considers the more complicated scenario where the new scalar couples to new colored 
particles, similar to [35j, as mentioned above. 

In addition to these renormalizable scenarios of NP effects on loop-induced SM Higgs phe- 
nomenology, a few papers have followed an effective field theory approach by constructing and 
constraining the size of dimension-six operators. In [7], the authors focused on the coefficients and 
constraints of operators for h — t- 77, jZ, and gg, while [12] extended the discussion to include 
h ^ ff as well. Importantly, both of these studies assume any New Physics contributions are 
heavy enough to be integrated out, thus there are no new particles in the low energy spectrum. 

In contrast to the previous literature, therefore, we discuss the general case using renormalizable 
interactions when both categories of NP effects are present. We isolate contributions with new 
colored scalars, new colored fermions, including Standard Model quark mixing, and new colored 
vectors, and we allow such effects to be modified by Higgs mixing. 

The paper is organized as follows. In Sec. |IT| we discuss general aspects of the Higgs portal 
relevant for our analysis of Standard Model Higgs production from gluon fusion. In Sec. Ill, we 



briefly review the leading order gg ^ h calculation for the Standard Model as well as the trivial 



extension of adding a fourth generation. In Sec. IV, we discuss gluon fusion in the presence of 
a new colored scalar. In Sec. IVl we present the analogous calculation for a general new colored 



fermion. Lastly, in Sec. VI, we discuss the interesting case of a new colored vector and its effects 
on gluon fusion in the context of a UV-complete, renormalizable model. Details of this calculation 



are presented in Appendix [A] We summarize and conclude in Sec. VII 



II. THE HIGGS PORTAL AND HIGGS MIXING 



In this section, we review the Higgs portal as a general framework for studying the connection 
between arbitrary New Physics models and Higgs physics, with a special emphasis on the resulting 
effect on gluon fusion. 

In the SM, the Higgs field is responsible for breaking SU {2)i x U{1)y gauge symmetry, resulting 
m masses for the and Z bosons as well as the chiral SM fermions. By virtue of being the 
only scalar field present in the SM, the Higgs also generates H^H, which is the lowest mass 
dimension operator possible in the SM that is both gauge and Lorentz invariant. Hence, arbitrary 



5 



NP operators can then be tacked on to H^H to give 

^hp ^ ^hpONpH'^ H . (1) 

Although Ojvp can be an arbitrarily high dimension operator, with an appropriate power sup- 
pression from a high scale A^rp, a generic Higgs portal term is only typically unsuppressed when 
Onp itself is dimension two and gauge and Lorentz invariant: hence, we take Onp ~ <I>^<I>. One 
exception is the case when a new scalar field is a pure SM and NP gauge singlet, but since we are 
focused on NP effects on gluon fusion, we will not discuss the gauge singlet case further. 

One class of NP effects on gluon fusion arises from new colored states that directly enter the 
gg ^ h loop diagram. The direct coupling of colored states to the Higgs via Eq. ([T| implies the 
mass of the new state is shifted after electroweak symmetry breaking (EWSB), and as this direct 
coupling is turned off, the NP effect vanishes. This class of effects is typified by models with new 
colored scalars, but a new fermion with Yukawa-like couplings to the SM Higgs boson also follows 
this scheme, albeit not via the Higgs portal. Although the case where the mass of the NP state 
arises primarily from the Higgs vev was discussed in j32j , in our more general framework the NP 
mass scale and the new couplings to the SM Higgs are independent. 

Since new particle masses do not have to arise from the SM Higgs vev, a second broad class of 
NP effects on gluon fusion emerges. Namely, if a new scalar field obtains a vev to spontaneously 
break a new gauge symmetry and if a Higgs portal term is present, this new scalar field will mix 
with the SM Higgs. In this class, NP effects coming from new colored states can infiltrate gluon 
fusion through the mixing induced from the Higgs portal even if these states do not directly couple 
to the SM Higgs. These effects will also exhibit the familiar non-decoupling features in SM gg ^ h 
loop calculations by chiral fermions or /i — )• 77 loop calculations by W bosons if the analogous NP 
states are present [I^IEZIEE]: however, this non-decoupling feature only applies to the new scalar 
field component of the scalar mass eigenstates. 

As mentioned in the introduction, we allow for both direct and Higgs mixing mediated categories 
of NP effects to be present simultaneously. These effects arise in many extended color sector models, 
and we consider isolated new colored scalars, fermions, and vectors in turn. For colored scalars, we 
couple them to the Higgs via the Higgs portal in Eq. Q, and hence they will exhibit an example 
of the direct category of NP effects with X^p as the direct coupling. For colored fermions, we 
consider two subcategories distinguished by the possibility of SM fermion mixing. If new fermions 
are introduced that mix with SM fermions, the usual SM calculation is modified to accommodate 
fermion mass eigenstates that do not typically couple with the SM Higgs with the usual Yukawa 
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strength. Without such fermion mixing, the SM calculation is unchanged and the new contribution 
arises from direct Yukawa couplings to the Higgs, the new scalar, or both. The decoupling behavior 
of new colored fermions are parametrized by fermion mixing angles and the possible scalar mixing 
angle. 

Perhaps the most interesting case is that of a massive colored vector boson. Here, in order to 
have a theory which is tree level unitary , it is natural to consider an extended color symmetry 
which is then spontaneously broken to SU{3)c gauge symmetry. Then the massive vectors corre- 
sponding to the broken generators form representations of the unbroken color symmetry. We are 
thus left with a renormalizable, unitary, spontaneously broken gauge theory |7Uj . 

We remark that another class of New Physics effects via the Higgs portal operator is possible. 
Broadly speaking, at the renormalizable level, where Onp ~ <1>^<I> in Eq. one class of Higgs 
portal effects is characterized by new colored scalars which do not obtain vevs. The second class is 
driven by new uncolored scalars that do obtain vevs from their scalar potential. Another possibility 
is colored scalars that do obtain vevs, but such color-breaking vacua are not viable phenomeno- 
logically. The last possibility consists of new uncolored scalars that do not obtain vevs from their 
scalar potential. Such a scalar does not enter the gg ^ h loop, but if X^p is large and positive, 
the resulting Higgs portal-induced shift in mass squared, —Xhpvf^ {vh is the Higgs vev) could drive 
the new scalar to acquire a vev. Hence, this last category of portal symmetry breaking models is 
unique because the Higgs portal coupling is a necessary ingredient for driving the new scalar to 
obtain a nonzero vev. Obviously, the roles of the new scalar and the Higgs scalar can be reversed, 
whereby the Higgs portal term allows a new scalar vev to drive the Higgs field to obtain a negative 
mass squared and hence trigger EWSB. We reserve a study of "Portal Symmetry Breaking" phe- 
nomenology for future work. Also, in the discussion above, we have delineated cases according to 
specific constraints on the Lagrangian parameters. A precise determination of these bounds would 
require an analysis of renormalization group evolution, which is beyond the scope of this work. 



A. New Physics Scalar Standard Model Higgs Mixing 

We briefly discuss the second class of NP effects from the Higgs portal described above, i.e. a 
new scalar and the SM Higgs both obtain vevs in Eq. ^ and mix. For simplicity, we only consider 
one new scalar, but our discussion is readily generalized to multiple scalars. We also assume $ 
transforms as a singlet under SU{2)i x f/(l)y, but that it is charged under a new local or global 
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symmetry in order to prevent "tadpole" terms. We let Oj^p ~ for a new scalar field giving 

£ D XhpH^H<^^^ ~ XhpVhvM ' (2) 
where we have suppressed representation indices and expanded the fields H ~ ~^J2^^ ^ ^^"^ 

<I> ~ —j={(t) + w^)- We assume the scalar potentials and V{H) are also present and Eq. (2) 

v2 IJ 

is the only Lagrangian term involving both <I> and H fields. The usual stability, triviality, and 
renormalizability constraints on the full scalar potential V{H) + V{^) — \hp\H\'^\^\'^ are assumed 



to be satisfied and will be imposed when we consider explicit models in Sees. |V] and VI Here, since 
^ obtains a vev, Eq. ^ leads to mixing via the mass matrix 

2 _ I < -\hpVhV4, \ 

"'-scalar " „ ' ^"^1 



I' 



-\hpVhV,l> m? 

where Vh and are calculated from minimizing the full potential V{H) + — A/jp|$p|-ffp 

and hence determine rrih and ra^. The functional dependence of m^^ and on their respective 
potential parameters can be fixed by solving the potentials V{H) and V{^) separately, and in 
the limit that A/jp — )• 0, the exact vevs Vh and recover their original, unperturbed values. This 
observation has important ramifications when calculating the exact Goldstone-Goldstone-scalar 



couplings needed for vector loop amplitudes in Feynman gauge, which are discussed in Subsec. |A 3 



We can readily diagonalize the symmetric mass matrix Eq. Q to obtain the mass eigenstates 

(4) 



si = hcosO — (psinO , 
S2 = h sin 9 + (j) cos 9 , 

with a Jacobi rotation mixing angle 9 defined by 



tan2g= -^y^l^ . (5) 



We will also need the inverse operations, 

h = si cos 9 + S2 sin 9 

4> = — si sin0 + §2 cosi 

The eigenvalues of Eq. Q are 

,2 _ 1 / 2 , _ ii. // _2 _^^2V ,A\2 „2„2 



(6) 



1 {rnl + ml)-\^ (-ml + m^) ' + AXl^y^ , (7) 



and 



m 



2 _ 1 1^2 I 2\ , ^ If ^2 I ^2^^ 



-i,yrnl + ml) + -d [-ml + ) + AXl^y^ , (8) 
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where we have taken rrig^ < without loss of generahty. As mentioned before and demonstrated 
in \27\ [301 135] , the mixing of the scalar states from the Higgs portal can significantly affect scalar 
production via gluon fusion. Moreover, the mixing is driven purely by the strength of X^p, which 
must be real but whose sign is not fixed. 



B. New Physics Effects on Production of si 2 

We can now readily disentangle the two categories of New Physics effects on gluon fusion. Now, 
because of h-cj) mixing via the Higgs portal in Eq. ([2]), we must calculate cross sections for gg — )• si 
and gg — S2 production instead of the gauge eigenstates h and (p. Since both h and cp can couple 
to new colored particles, contributions to gg — )• si^2 can manifest themselves through both the h 
and (p components of ,31^2; leading to suppression or enhancement of the production rate relative to 
the SM. This also implies that partial decay widths are affected, whereas in hidden sector models, 
such widths are unaltered apart from a universal cos^ 9 suppression coming from Higgs mixing. 

From the discussion above, we can decompose the production amplitude of si via gluon fusion 
in terms of the gauge eigenstate h and (p production amplitudes as. 



Migg ^ si) = ce [M{gg /i)]L^=^^^ se [M{gg ^ Mr 
M{gg S2) = se [M{gg /i)]L,^=^^^ +ce [M{gg - 



(9) 



where cq = cos 9, sq = sm.9 are defined by Eq. ([5]). In the discussion below, we presume the matrix 
elements are evaluated at the appropriate scalar mass and will drop the notation above. Hence, 
given the linear combination dictated by Eq. ([9]), we are now free to isolate the contributions to 
gg ^ h and gg (p. 

We are particularly interested in identifying, at the amplitude level, the mechanisms respon- 
sible for modifying gluon fusion and whether and how they can decouple. A completely general 
expression for all possible NP effects along these lines is cumbersome, so instead we write 



M{gg si) = ce 



se 



M{gg > h) + M{gg > h) + M{gg — h) 

scalars fermions vectors 

M{gg — - — > (p) + M{gg — > 4>) + M{gg 



(10) 



scalars fermions vectors 

and treat each category of loop particles separately^ Each of these categories can be further 
subdivided into particles that couple solely to /i, solely to (p, or simultaneously to both. In the 



^ For the vector loop calculation, we implicitly assume a unitary gauge calculation where only vectors appear in the 
loop. If working in Feynman gauge, the associated Goldstone and ghosts would also be part of the vector category. 
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scalar case, for example, we can write 



M{gg 



scalars 



Si} 



ce 



se 



Y,M{gg-^h) + Y.M{gg-^h) 



EMigg 



+ EM{gg-^ 

k 



(11) 



where the scalars rji, rjj, rjk couple only to h, both to h and (p, and only to (j), respectively. We can 
now make definitive statements about the decoupling behavior of the scalars r]i, rjj and rjk- If the 
masses of rji [r]/.) arise solely from the vev {V(f>), then these scalars will exhibit non-decoupling 
from h {(j)) as their masses are taken very large: if instead their masses include sources besides 
Vh or V(fi, then decoupling will occur as the mass scale of these new sources is taken large. The 
behavior of the r]j states are a straightforward combination of the previous arguments. 
For fermions, we write 



M{gg 



fermion. 



■> Si) 



Ce 



se 



^Migg-^h) + ZMigg-^h) 



EM{gg^ 

j i>j 



+ EM{gg^ 

k 



(12) 



To be more illustrative, we can take some familiar examples to demonstrate the flexibility 
of Eq. (12). In the case with Higgs mixing but without new fermions ij^j or ^fc, then ipi con- 



sists of the SM quarks and we get a universal cg suppression of the matrix element. If instead we 
only add a new vector-like top partner to the SM, then = 1, = 0, and ipi includes the first 
five SM quarks and the two fermion mass eigenstates resulting from top mixing while the ifjj and 
ipk sums are absent. Finally, if Higgs mixing is present and new colored fermions are added that 
couple both to h and (f) but do not mix with the SM fermions, then ipi will run over the SM quarks 
and ifjj will run over the NP colored fermions. 

Lastly, we can introduce massive colored vectors. We will only consider the case where these 
vectors couple to (j), giving the relatively simple expression 



M{gg 



vectors 



> Si) 



-Se 



EMigg-^ 

k 



(13) 



emphasizing that this contribution to the gluon fusion rate for si production relies on the Higgs 
portal, since the SM Higgs is assumed to play no role in breaking the extended color gauge sym- 
metry. 

After the above discussion, we present a parametric understanding of how production and decays 
of si 2 are affected by direct coupling and h-(j) mixing. As we have seen, performing a completely 
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general analysis would be overly cumbersome, and so we will make a few mild assumptions to make 
the analysis more intuitive and tractable. Throughout the discussion, we assume a narrow width 
approximation, allowing us to factorize production and decay processes. 

We define the overall leading order enhancement or suppression factor of si production relative 
to SM Higgs production via gluon fusion as 



c^igg si) \Migg si)|^ \ceM{gg ^ h) - seM{gg 



^ SM ' 



M(gg > h) 

SM ' 



Migg )■ h) 

^ SM ' 



Cq \ Zggh tgZgg^l 



using Eq. ^ and with tg = tan 9. The complex amplitude ratios are given by 



Z. 



ggh 



Mjgg^h) 

Migg > h) 

^ SM ' 



-'gg<t> 



M{gg^4>) 
M{gg > h) ' 

^ SM ' 



(14) 



(15) 



and will simplify significantly for any given NP model under consideration, as we will demonstrate 
in Sees. IV, Vjand VI We see that both egg > 1 (signaling enhancement) and egg < 1 (signaling 
suppression) are possible with New Physics and changing the sign of A/ip. In the limit that 6 = 0, 
the only effect on gluon fusion arises from the inclusion of new colored states that directly couple 
to the SM Higgs, which was a main focus of \19\ [20l I32j . In the case where Higgs mixing is the 
only new effect, then Zggh = 1 and Zgg^f, = 0, and we have the simple expression egg = Cg, as noted 
in [27]. 

We remark that complete suppression of gluon fusion does not correspond to vanishing LHC 
production for the si state. This is because the subdominant modes of vector boson fusion, vector 
boson association, and tih production comprise 12.5% of the total cross section for a SM Higgs 



mass at 125 GeV [71J. Moreover, even if the leading order cancellation in Eq. (14) is exact, we 
expect NLO corrections, which can be as large as 20% in the case of colored stops [8j, to make the 
cancellation imperfect. 



C. New Physics Effects on Decays of si 2 

We now extend our discussion to include NP effects on decay widths for our scalar state si, 
which we take to be dominantly SM Higgs-like. We will not detail all of the (practically infinite!) 
possible final states for si, but will instead focus on the WW, ZZ, 77, bb and r+r^ decay channels. 
For the WW final state, we write 

M{si WW) = ceM{h WW) - seM{(t) WW) w ceM{h WW) , (16) 
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and thus 



B(h >WW) "Ts, 

SM 

where we have assumed the tree-level coupling of hWW dominates the (typically loop-induced) 
coupling of (pWW, and F/j and F^^ are the total width of the purely SM Higgs and the mass 
eigenstate si, respectively. Under the same assumption that hZZ dominates the (f)ZZ coupling, 



the same result in Eq. (17) also applies to the ZZ final state, and so branching ratios of si to WW 
or ZZ diboson states are typically suppressed in Higgs mixing models. 

For the diphoton final state, we can adapt our gluon fusion discussion, replacing colored particles 
with electromagnetically charged particles. Following the guide of Eq. ([£]), this gives 

M{s, ^ 77) = ce [Mih ^ 77)]L,=^.^ - se [MicP ^ 77)]L,=^.^ • (18) 
Unlike the WW or ZZ decay modes, the /i — t- 77 decay is induced at loop level in the SM and 



new contributions can easily cancel against or add to the SM contributions. Using Eq. (18), we 
can write the relative branching ratio as. 



where e^^ is analogous to egg in Eq. (14) and Z^^^ and Z^f,^^ are defined similarly. 



The relative rate for gg ^ si ^ 77 versus ^(g — )• /i — )• 77 is now given by 



T 

i Sl 



In many models, though, the various inputs for Eq. (20) reduce to simple expressions. For example, 
in Higgs mixing scenarios where (f) only couples to hidden sector particles, we obtain Zh'yy = 1, 
^(^77 = 0, and so 

which agrees with the universal Cq suppression noted in [27]. Another simple limiting case arises 
if we take 6 = and introduce new charged particles in the 77 loop coupling to the Higgs. In this 
case, h = Sl and we can write 



SO that only the direct NP effects contribute. 
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Finally, we can calculate the si branching ratio to bb or r+r . If Higgs mixing is present, if 
does not appreciably couple to the SM fermions, and if the SM fermions are not mixed with 



NP fermions, then the same results from Eq. (17) apply, substituting // for WW. A completely 



general expression, however, because of the possible presence of all of these effects, is unwieldy. As 
an explicit case, for the bb final state, if we allow for h-(j) mixing and introduce a coupling between 
(j) and bb, we obtain 

SM 

where = ^ bb)/M.{h > bb). We can see that interference effects from Z^f^, although 

SKI 

diluted by tg, can lead to an overall increase in the branching fraction of si — t- bb. 

In summary, we have disentangled the effects from Higgs portal-induced mixing of h and (j) 
from NP effects caused by direct coupling to /i, (/>, or both. For gluon fusion, we have explicitly 



identified the decoupling behavior of New Physics states in Eqs. (11), (12) and (13). If we assume 
NP couplings to be small or negligible, then the resulting si branching ratio has a universal 
Cg suppression and a universal total width ratio suppression. On the other hand, interference 
effects resulting from couplings to h and/or (f) lead to a myriad of effects and possibilities for both 
suppression and enhancement of relative rates can be achieved. 

We note that all of these expressions can readily be adapted for S2 decay with an appropriate 
C0 — )• se, —S0 —7- cg exchange and m^^ — )• m^j. In addition, if > ^mg^, there is the additional 
decay mode S2 — )• sisi, as emphasized in [30]. Also, if any of the new states are lighter than 
or 772^2/2) then additional non-standard decay modes open up. This effect is manifest in the above 
expressions through the ratio of total widths F/j/Fs^. 



III. THE gg^h PROCESS IN SM 

Here we briefiy review the leading order Standard Model calculation for Higgs production via 
gluon fusion. As shown in Fig. [T| gluon fusion arises in the SM via quark loops, with the dominant 
contribution coming from the top quark with its large Yukawa coupling. We again highlight the 
fact that since neither the W or Z boson couplings are probed in this production mode, large 
effects can be present in this loop process that strongly change Higgs production but do not affect 
EWSB. 
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FIG. 1. The Standard Model contribution to gg — h. 



The total Standard Model amplitude is 



f f 



TT 



Vh V ^ 



(24) 



where / runs over the SM quarks, C{rf) is the Casimir invariant {C{rf) = 1/2 for SM quarks), 



a and d are color indices, pi ■ p2 
well-known loop function 



ml 



for an on-shell Higgs, rj = m|/(4mj) and Fp{t) is the 



Ff{t) 



(r + (r-l)Z(r)) , 



(25) 



using 



Z(t) 



log 



arcsm ^JT 
1 + Vl -r-i 



— 1 



T < 1 
- 2 



r > 1 . 



(26) 



1 - y/l- T- 

Because the SM quarks obtain their mass purely from EWSB, they do not decouple even as we 
take the Umit r — >• (equivalent to mj — )■ oo). In the case of the SM4, this sum would include t' 
and b' contributions. In the limit that the SM Higgs is well below the threshold for t, t' , and b' to 
propagate on-shell in Fig.[T| we obtain the usual factor of 3 enhancement in the SM4 matrix element 
for gg — )■ h, which gives, at leading order, a factor of 9 enhancement for gluon fusion production 
in SM4 relative to SM3 [56j. Recent literature that has attempted to resolve the quandary of 
a sequential fourth generation of fermions with the lack of enhancement in gluon fusion include 
Refs. inUETHBO]- 

We can anticipate, in the presence of new additions to gluon fusion coming from New Physics, 
that the main structure of Eq. ( |24[ ) will not change apart from possible differences in the scalar 
vertex, C(r), and the loop function F{t). In particular, the P1P2 — Pi ■ P2g^'^ structure of the 
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matrix element is assured by 5C/(3)c gauge invariance (or equivalently, the Ward identity). This 
is analogous to the situation in the /i — >• 77 calculation, where electromagnetic gauge invariance 
requires the same momentum structure |68j . 

IV. NEW COLORED SCALAR 

In this section, we isolate and calculate the effect of a colored complex scalar S propagating in 
the gg ^ h loop. We use the Higgs portal in Eq. ([2| to couple S to the SM Higgs, and we write 
a (positive) tree-level mass squared for S such that SU{3)c remains unbroken and Higgs mixing 
is absent. Depending on the sign and strength of Xhp, we can achieve significant suppression or 
enhancement of gluon fusion as a result of the interference between the SM fermions and the colored 
scalar. 

The Lagrangian involving S is 

Cs = \D^S\'^ - mlS^S - k\S^S\'^ + XhpS^SH^H , (27) 

where color indices have been suppressed and we assume mg > and k > to ensure stability. 
As discussed in Sec. [ll| Xhp must be real: for positive (negative) A^p, we will get destructive 
(constructive) interference with the SM loop calculation, in agreement with |20^ [35] (note we have 
a different sign convention for Xhp)- After EWSB, the physical scalar mass is 

"T-l = "ig - Xhpvl , (28) 

which imposes the constraint that rn-Q > Xhpvf^ to avoid portal symmetry breaking of SU{3)c- 

The two diagrams to calculate are shown in Fig. [2] Since S is complex, the matrix element 
for Fig. [2]A needs to be multiplied by 2 to account for the charge conjugate diagram: if S were 
real, no factor of 2 is used and instead the matrix element for Fig. [2p must include a symmetry 
factor of (1/2). 

The total amplitude corresponding to the diagrams in Fig. [2] for a complex scalar field propa- 
gating in the loop is 

'•^5" = • (?) - ^S^Fsirs) . (29) 

where C{rs) is the SU{3)c Casimir invariant for S, ts = m|/(4m|) and the loop function Fs is 
defined to be 

Fs{t) = - t-^Z{t) , (30) 



15 



s I 



(A) 




(B) 



FIG. 2. Feynman diagrams for scalar foop contributions to gg — > h. For a complex scalar one must also 
include the charge conjugate equivalent of diagram (A). 



with Z{t) defined in Eq. (26). Note that for fixed X^p the amplitude decouples as ms 
of the tree- level mass squared tjiq. 



oo because 



Now, the summed amplitude for M{gg 



SM+S 



> h) is 



m 



/ 



Mr/. 



(31) 



If ms, rrit > mfi/2, then Fs is strictly real and negative and Fp is strictly real and positive, which 
implies that for A^p > {Xhp < 0) the interference between the colored scalar amplitude and the 
SM is destructive (constructive). 

Since the phase space integration needed to calculate the si cross section is identical to the SM 



Higgs case, we can write the ratio from Eq. ( 14 ) as 



f V Vh 



C{rs)Xh pVh 
4m c 



"aaiSM+s 



(32) 



We consider the addition of a real color octet scalar (C(rs) = 3, symmetry factor of 1/2) and 
a complex color triplet scalar (C(rs) = 1/2) and plot egg in Fig. |3] as a function of ms for some 
representative choices of parameters mh and Xhp. For the SM calculation, we sum over bottom 
and top quark contributions, using m^ = 4.20 GeV and mt = 172.5 GeV. We adopt the results 
of |72j to draw vertical exclusion bands on light color octet scalars from ATLAS [T^ and CMS |74j 
in dijet pair resonance searches. The gap in sensitivity from 200 GeV to 320 GeV between the 
34 pb^^ ATLAS search and the 2.2 fb^^ CMS search is a result of the increased multijet trigger 
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150 200 250 300 

ms(GeV) 



350 



400 



SM + Color Octet Scalar (m/, = 500 GeV) 
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FIG. 3. Relative rate egg in Eq. (32), showing the effect of the inclusion of a real color octet scalar (top 



row) or complex color triplet scalar (bottom row), for m/i — 125 GeV (left panels) or m/,, = 500 GeV (right 
panels). At the center of each plot, from top to bottom, the solid lines correspond to Xhp = —0.2, —0.05, 
0.05, 0.2. The left (right) gray bands in the octet scalar plots come from the ATLAS (CMS) search for 
pair produced dijet resonances. For the triplet case, the CMS bound still applies but the ATLAS bound is 
unconstraining after rescaling cross sections. 



threshold to handle more difficult run conditions. In particular, for the CMS study, turn-on effects 
of the QCD multijet trigger made the background modeling unreliable below 320 GeV. 

We see that both enhancement and suppression are possible, typically delineated by the choice 
of the sign of Xhp- This is manifest in the region where ms > rnh/2 where Xhp > corresponds 
to a suppression and Xhp < corresponds to an enhancement. In the region where ms < mh/2, 
we see enhancement and suppression for both signs of Xhp since the scalars can go on-shell in 
the loop, leading to an additional imaginary contribution to the scalar amplitude. The resulting 
interference is complicated by our inclusion of the bottom quark and its imaginary contribution, 
so the overall magnitude has competing cancellations among real and imaginary amplitude pieces. 
We note that Fig. [3] shows the expected decoupling of S as ms grows. We also remark that for 
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negative values of Xhp, our results are consistent with [20], where the finite difference in our results 
is a result of our inclusion of the bottom quark. Lastly, with regards to the ATLAS dijet pair 
search, we note that the complex triplet scalar is 1 /9 the production cross section of the real octet 
scalar, if their masses are equal, rendering the search insensitive to complex triplet scalars. 



V. NEW COLORED FERMION 



Adding new colored fermions to the Standard Model can greatly affect gluon fusion SM Higgs 
production in a number of unique ways. On one hand, new sequential generations of chiral fermions 
will add constructively with the SM fermion loops and, at leading order, scale the top quark loop 



by a multiplicative factor, as discussed in Sec. Ill On the other hand, a new vector-like colored 



fermion that does not mix with SM fermions has no effect on gluon fusion. In general, a new 
colored fermion mass eigenstate comprised of chiral and vector-like components will enhance the 
SM Higgs gluon fusion rate according to the chiral projection of the mass eigenstate. 

Because we also allow for Higgs portal-induced scalar mixing, though, the general situation can 
lead to either an overall suppression or enhancement of the gluon fusion rate. A model demon- 
strating the myriad of competing effects is straightforward to construct but only illuminating in 
its limiting cases. Hence, we will initially consider only mixing between a NP fermion and a SM 
fermion, neglecting the Higgs portal and Higgs mixing. 

We add new vector-like top partners [THl EB] , given by 

XL,R ~ (3, 1)2/3 • (33) 

This leads to the following mass terms, 

C D -ytHQitR - VlHQlXr - MxLXR + h.c. , (34) 

where M is a free parameter and induces mixing between the SM top quark and x- the (t, x) 
gauge basis, we have mass M and interaction Nh matrices given by 



^h={ , (35) 





with = — and Mf = — Note the entry in M can always be ensured since it corre- 
V 2 v 2 

sponds to the M'xhtR operator, which can be trivially rotated away since XR and have the 
same quantum numbers. The mass basis rotation matrices are defined in the usual way from 
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Md 



and L{MM'^)D 



Md 



After rotating M and Nh on the left (right) by a 



left-handed (right-handed) rotation matrix, we obtain 



/: D -f ( Mfl + —Vh ]PRt+ h.c. , 

Vh 



(36) 



where t = (ti,t2) and Md = LMW = diag(mt^, mt2), Vh = LNhR^. The coupling matrix Vh 
dictates the couplings of the top sector to the SM Higgs and, in principle, can have off diagonal 
entries; however, only the diagonal entries contribute to gg — >• h, because the L and R rotations 
leave the gauge interactions diagonal in the mass basis. 



Pi — ^ n,a 
9 ^751515151515151T>- ^ 




FIG. 4. Exotic fermion contribution in the mass eigenbasis. 



In this top partner scenario, each mass eigenstate gives its own contribution to the diagrams 
in Fig. [4| Since SU{3)c gauge invariance guarantees these two contributions differ only in their 
coupling to the Higgs via the element Vh, we can take the SM result for gg ^ h through 
a top quark and insert the appropriate element of Vh in place of the usual Yukawa coupling. 
This approach also encompasses more complicated mixing scenarios whereby the matrix element 
will exhibit different combinations of mixing angles and couplings for the various V/i entries as a 
replacement for the hff vertex in the gg ^ h amplitude. Since we are focused on exotic fermion 
effects on gg — )■ h, we take the (yh)ij entry to be a free parameter, which can be readily calculated 
in any complete model. 

The amplitudes involving exotic fermions in the loop are 

iMf = iV m ( ^ ( ^) (pW - 4,"^) FHrp.) , (37) 



Vh J V 2 



where the repeated indices on {Vh)ii are not summed, Ff{t) is given by Eq. (25 ), Tp^ = m^/ (4m|..), 
and Fi £ {ti,t2}- We see that the amplitude decouples as mp- — )• oo, unless mp^ and (yh)ii are 
generated by a common scale such as the Higgs vev. These direct new physics contributions 
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will alter gg ^ h even in the absence of Higgs mixing. Generally these contributions will add 
constructively if {Vh)ii > 0. 

Now, we augment the previous discussion to include Higgs mixing between h and a new scalar 



. We replace the vector- like mass M in Eq. (34) by a Yukawa term which generates the desired 



mass term and a new interaction term involving (f) after (p obtains a vev, giving 



y<j><PXLXR =^ M{1 H )XLXR 



(38) 



where M 



We now have a second interaction matrix V^, which is added to Eq. (36) and 



defined analogously to Vh, where 




(39) 



As discussed in Subsec. JIB the ti, t2 loops will give an enhancement factor for gluon fusion 
production given by 



''99\SM+XL,R ~ 

'C{rf) 



E 

/, no t 



Vh 



E ( 



Vh 



(40) 



Note that as mt^^t^ — ?■ oo, ti will decouple from the (j) component of si but not the h component, 
and vice versa for t2. 

We recognize that these new top partners will induce shifts in the electroweak oblique parameters 
S and T, but a full analysis of the allowed top partner parameter space is beyond the scope of this 
study. Therefore, we adapt the results from Ref. jl9], which studied top partner effects on Higgs 
production and included the constraints from the S and T oblique corrections. We set > m|, 
and for the vevs we are considering, —XhpVhV^ is a small perturbation to the diagonal mass terms 
in Eq. ([s]), allowing us to approximate the si contribution to S and T by the Higgs contribution 
considered in [19j for equal masses. We can thus illustrate our main point, suppression of gluon 
fusion, in this phenomenologically viable top partner scenario. In Fig. [5] we plot contours of €gg as 
a function of the left-handed fermion mixing angle and the heavy fermion mass eigenstate, for 
representative values of Xhp- The shaded bands correspond to regions consistent with the oblique 
parameters at the 68 % and 95 % C.L., taken from [19]. 
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Separately, we can also consider new colored fermions which do not mix with the SM quarks. 
As a final example, we consider new electroweak singlet fermions in the adjoint and fundamental 
representations of SU (3)c which obtain mass from the new Yukawa term of Eq. ( |38[ ) (with x — )• -0). 
In particular, these fermions do not couple to h, and hence their effects will be suppressed by the 
scalar mixing angle in Eq. ([s]). The relative rate is now 



e(^.h., 



te ('^F^(r^) 



( Cjrf) 



(41) 



We can see that ip does not decouple from the gg — )• si amplitude as its mass is taken very 
large because Fp{t^) asymptotes to a finite value. We show egg in Fig. [6]for two choices of color 
representations. We see from Fig. [ojthat the octet fermion (which includes a 1/2 to account for 
lack of conjugate diagram for a real fermion) produces larger suppression or enhancement than the 
triplet fermion for identical Xhp values, as expected from the difference in their respective C{r^). 
Collider constraints on these new fermions will require model dependent assumptions about their 
decay channels, and since our focus is on the model independent effects on gluon fusion, we do not 
consider such constraints here. 



SM + top partners SM + top partners 
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■ 500 1000 1500 2000 " 500 1000 1500 2000 

m,,(GeV) mfc(GeV) 

FIG. 5. Contours of the relative rate of si production as a function of rrit^ and the left-handed mixing angle 
in the top partner scenario for = 800 GeV, = 500 GeV, and Xhp — —1 (left) and 0.75 (right). We 
adapt the analysis and results of Ref. [19] to show shaded contours which are consistent with the oblique 
parameters S and T at the 68% and 95% C.L.. 
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SM + New Dirac Fermion (8,1)o SM + New Dirac Fermion (3,1)o 
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FIG. 6. Relative rate of Si production with Higgs mixing and a new color octet fermion (left) and a new 
color triplet fermion (right). We choose — 500 GeV, m^^ = 500 GeV and m^j = 700 GeV. From top to 
bottom, the solid lines in each plot correspond to A^.^ — 0.015, 0.005, —0.005, —0.015. 



VI. NEW COLORED VECTOR 



The last type of New Physics contribution to gluon fusion we will consider is the addtion of a 
new massive colored vector. In a renormalizable theory, the massive vector must arise from a spon- 
taneously broken gauge theory, which necessarily entails the addition of a new scalar that acquires 
a vev and can mix with the SM Higgs via Eq. ([2|. Correspondingly, we will consider an extended 
gauge symmetry SU{3)i x 5C/(3)2, known as the renormalizable coloron model (ReCoM) \63\ 164]. 
In this model, the complex scalar field <I> transforms as (3, 3) and obtains a diagonal vev (when 
written as a matrix- valued field), which breaks SU {3)i x SU {3)2 to the diagonal subgroup, which is 
identified with the SM SU (3)c symmetry. The field has 18 degrees of freedom: 8 are "eaten" by 
the broken gauge generators to make the massive color vector G' known as the coloron, 8 become a 
real scalar SU{3)c octet Gh, and the remaining 2 are the real scalar 0/j and pseudoscalar color 
singlet fields. Hence, in this construction and a consequence of the Higgs portal, the addition of a 
massive color vector G' concomitantly includes a new scalar octet and two new scalar singlets, of 
which Gh necessarily affects gluon fusion and (pji mixes with the SM Higgs. 



A. The Renormalizable Coloron Model 



We analyze the total scalar potential including the SM, the ReCoM, and the Higgs portal. Our 
analysis mirrors that found in [HI] , except our addition of the Higgs portal operator creates a link 
between the two scalar potentials V{H) and V{^) and hence the two vevs must be solved for 
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simultaneously. The full scalar potential is 



Vtot = y{^) + V{H) + Vhp 



A 



Tr($t$) _^^(det$+ h.c.) + ^ [Tr($^>t)]^ + ^ Tr(^>$t$$t) 



(42) 



- \hp\H\'' Tr($t$) , 

where we assume > without loss of generality. We require ni\ > and > such that $ 
and H will acquire vevs. We also require 3A$ + k<j> > and A/^ > so each individual potential is 
bounded from below: we neglect renormalization group effects when discussing bounds on potential 
parameters. 

It is straightforward to find the vevs for $ and H by decoupling the two equation system. We 
find, in analogy with 




2 , Xhprrijj 
Xh 



2(3A$ + K$) 



-I3 . 



(43) 



If A/ip is too large, then it can generate a new ground state at large field values of h and 0. The 

2 

resulting upper bound on X^p is A^^ < -Ah(3A$ + «;$), which we satisfy by requiring > 0. 



Given Eq. ( 43 ) , the Higgs vev is most easily written as 

/ 

1 / \ 1 




1 

71 







V 



m 
'Xh 



(44) 



2Xh I 



and Vh is fixed to be 246 GeV. 

Expanding $ around its vev, we get 

1 



CD = — (^^ + 0^ + ic^i) I3 + (G\ + iG'^c) 



(45) 



where and are S'?7(3)c singlets and Gh and Gg are SU{?>)c octets [M]- The Gq comprise 
the Goldstone bosons which give mass to the coloron, G' . The Higgs is decomposed in the usual 
way. 



H 



1 



^ \ vh + h + iGo 

where Go and G^ are the Goldstone bosons eaten by the electroweak gauge bosons. 



(46) 
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After the spontaneous symmetry breaking SU{3)i x SU{3)2 — )■ SU{3)c and EWSB, mixing is 
induced between the singlets (pn and h. This leads to the mass squared matrix in the {h, 
interaction basis given in Eq. (|3|) but with m| — t- m?, and 



ml = 2XHvl, m^^ = |(3A<, + K<,)-^, (47) 

where Vh and depend on Xhp. The assumption of u ^ > and the conditions < m\ < 
imply 

< y^C^-^* + i^<s>)v^ ■ (48) 



(49) 



By our assumptions, the right hand side of Eq. (48) is positive definite and thus bounds fi^ from 
above. Our analysis follows exactly the same steps as Subsec. II A giving the following results: 

tan2e = Z2AW 
TTij, — mi 

51 = hcosO — (j)Rsm6 , 

52 = h sin + (pfi sin 9 . 

For the physical masses of si and S2, we obtain 

For the physical masses of the remaining scalars in the spectrum we find, 

'^li = ^J^f^'S'Vcj, , rnl;^ = ^{2ml^ + k$u|) , (51) 

which agrees with |64| in the limit Xhp — >• 0. 

The constraint m^^ > is satisfied since we assumed fj,^ > and > 0. Requiring rri^^ > 



2 

ml 
si 2' 



,2 



(50) 



implies m^^ > which augments the previous condition Eq. (48) to give 



V6 ^ ^ V 3 
In order for a valid range of /j,^ to exist, we thus require 



K^V^ / 2 

< /i$ < A / - (3A<i, + K$) (52) 



2A$ + > . (53) 



Our subsequent analysis ensures these constraints are satisfied. 
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B. Phenomenology 



The diagrams for the colored vector loop in unitary gauge are shown in Fig. [7| Although the 
ReCoM model includes the Gh scalar octet contribution, we have isolated colored scalar amplitudes 



in Sec. IV, and so we focus here on the colored vector contribution. 




(A) 




(B) 



FIG. 7. Feynman diagrams for vector loop contributions to gg — > (pR in the unitary gauge. 



The full amplitude for a real vector field propagating in the loop is 



where 



\ms„=ms 



<t>R — "-''l 



(54) 



Fv{t) = - (r-i(3 + 2r) + 3r-2(_i + 2t)Z{t)) . 



(55) 



A full derivation of this amplitude in both unitary and Feynman gauge is presented in Appendix]^ 
The resulting enhancement factor is 



"gglSM+ReCoM ~ 



( Cjrf) 



(56) 



where xg^j/i 



a + i™'^/ ~ '^^Gh] /^<^ evaluated at m^^ = m^^ in Eq. (56) is the Gh 



coupling to (pji and we have included both Gh (with an explicit symmetry factor of 1/2) and 
G' contributions. We note that the real colored vector loop function Fqi is numerically about 
a factor of 5 larger and of the opposite sign than the usual SM loop function Fp- The scalar 
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loop function Fs{t) is also of the opposite sign and roughly a third of Fp: the loop functions 
are plotted in Fig. [Sj We comment that, as a result of the large loop function for the colored 
vector, moderate values of Xhp can have large effects on gluon fusion production. This provides a 
straightforward construction, for example, to counteract the enhancement from a fourth generation 
of chiral fermions in the Standard Model. If such a large cancellation oi gg ^ h amplitudes was 
present, however, we expect di-Higgs production via gg — t- hh to be correspondingly enhanced if 
the gg — )• hh triangle and bubble amplitudes are negligible: this is because the individual signs of 
direct Higgs couplings that lead to suppression become squared in the gg — hh box amplitude. In 
this case, the di-Higgs gluon fusion production channel may be a promising discovery mode, and a 
more careful study is certainly warranted. 




FIG. 8. Loop functions — Fp (top), Fl — Fs (middle), and F^ = Fy (bottom) for fermion, scalar, and 
vector particles of mass M, respectively. The loop functions develop an imaginary part for r > 1, which 
corresponds to the particles in the loop going on-shcU. 

We first present results in Fig. [9] for Cgg with the sole addition of the colored vector for i;^ = 200 
GeV and 3 TeV. The = 200 GeV choice is interesting because bounds on colorons in the 
low mass region coming from pairs of dijet resonances from ATLAS [73j and CMS [74j are not 



currently sensitive in the 200 to 320 GeV range, as discussed in Sec. IV There are numerous recent 
studies which have done fits of the Higgs couplings to the data, including the Higgs-gluon effective 
coupling [TTHSU] . For v^p = 200 GeV, a coloron mass of 250 GeV and mixing angle of sg ~ —0.08 
one can reproduce the best fit point {cg ~ 0.5) for the scalar-gluon effective coupling found in [78] 
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and as can be seen from Fig. [9] one can also easily reproduce values which give rates close to the 
SM value should the effective coupling migrate towards the SM prediction with more data. We 
also consider = 3 TeV, which is the scale probed in dijet resonance searches using 4 fb^^ of 
8 TeV LHC data at CMS [81j- The observed limit from this search on the coloron mass, mQ/, 
is 3.28 TeV. The bounds are indicated by gray vertical bands as described in the caption. The 
latest dijet resonance search done by the ATLAS collaboration [82] with 5.8 fh'^ of 8 TeV LHC 
data does not report an observed limit on the coloron mass but we expect the limit to be just 
a little higher because of the increased luminosity. An additional constraint on the coloron mass 
arises from the constraints imposed by ReCoM, whereby perturbativity restrictions on the original 
SU{3)i X SU{3)2 gauge couplings give an upper limit and requirements on generating the correct 
SU{3)c coupling give a lower limit. In deriving these bounds, which are given by dashed vertical 
lines in Fig. |9j we have neglected renormalization group running of as- 

In general, the ReCoM model includes contributions from the color vector G' and the scalar octet 



Gff. We can see from Eq. (56) that the contribution from Gh coming through the h component 
of si always adds constructively with the G' contribution. Whether the contribution from Gh 
entering through (p^ also adds constructively with the G' contribution depends on the sign of xg^, 
which in turn depends on the mass hierarchy between the various scalars. We present the effect 
arising from only the color vector in the top row and from both new colored states in the bottom 
row of Fig. [9j For these plots, we have set m^i = 125 GeV, = 225 GeV, m,^^ = 160 GeV, and 
= 140 GeVg 

We remark that the flat behavior of €gg in each plot arises because for mc > mg-^, the loop 
function dependence of niQi asymptotes quickly. This reflects the fact that as niQ/ is taken large, its 
effects (which enter only through the (pR component of si) do not decouple from the si production 
amplitude, which is reminiscent of the non-decoupling of W bosons from the SM Higgs in /i — )• 77. 
In addition, for the scalar octet Gh, which couples to both h and <j)Ft components, we find that as 
mcjj is taken large, its effects do not decouple from the (j)R component but do decouple from h, 
see Sec. IV Finally, we note the small reduction of Cgg in ReCoM is a result of the Gh contribution 
slightly cancelling the G' contribution, given our chosen parameter point for which < 0) 

and we see that the overall effect is dominated by the coloron contribution, as expected from the 
magnitudes of the loop functions shown in Fig. [8) 



^ The ATLAS search for dijet pairs has an upward fluctuation above their expected limit in the 140 GeV range, 
which leaves the window of a scalar octet open for this mass point. 
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FIG. 9. Relative rate of si production showing the effect of adding only the coloron (top row) and the 
effect of adding the coloron and Gh, with toq^ ~ 140 GeV (bottom row). From top to bottom, the solid 
lines correspond to Xhp — —0.0015, —0.001, 0.001, 0.0015. In the left panels, the gray bands correspond 
to the pair-produced dijet bounds from ATLAS (left) and CMS (right). In the right panels, the left gray 
band corresponds to the CMS exclusion on coloron production in dijet resonance searches. The left vertical 
dotted line corresponds to the minimum mc mass allowed given the specified choice of v^, and the right 
vertical dotted line corresponds to the perturbativity constraint. 



VII. DISCUSSION AND CONCLUSION 



We have seen that new colored particles can have significant effects on gluon fusion production. 
Also, the inclusion of Higgs portal-induced scalar mixing readily leads to new possibilities for 
suppressing or enhancing the gluon fusion rate. We have isolated contributions arising from new 
colored scalars, new colored fermions, including quark mixing, and new colored vectors. With the 



amplitudes in Eq. (29), Eq. (37), and Eq. (54), we can immediately calculate the interference effects 
present in general new physics models. We have demonstrated these effects can easily run from 
0{1%) to 0(10%), and both suppression and enhancement of gluon fusion can occur. Moreover, 
such large deviations are possible by colored states at mass scales that can be directly probed at 
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the LHC. In particular, when the effects on gluon fusion scale with the Higgs portal coupling A/ip, 
the dearth of restrictions on X^p from the underlying theory can lead to very large modifications. If 
many competing effects are present in the gg ^ h amplitude, we expect that di- Higgs production 
will be correspondingly altered and a promising discovery mode, but we leave a more careful study 
for future work. 

Since gluon fusion production does not directly probe the Higgs mechanism or the phenomenon 
of electroweak symmetry breaking, Higgs identification studies should generally allow for mixing 
with a separate scalar state as well as competing effects from multiple new colored states running 
in the loop. Our general framework and analysis can be easily mapped onto the various recent 
studies which attempt to fit the data with non-SM Higgs couplings to two gluons. In particular 
we have shown that the various new physics effects can conspire to give rates very close to the SM 
expectation as well as easily accounting for any slight deviations as suggested by recent fits of the 
scalar effective couplings [77ti80] . 

ACKNOWLEDGEMENTS 

The authors are grateful to Bill Bardeen, Marcela Carena, Andre de Gouvea, Bogdan Dobrescu, 
Jennifer Kile, Ian Low, Adam Martin, Reinard Primulando, Pedro Schwaller, Daniel Stolarski, and 
Tim Tait for useful discussions. Fermilab is operated by Fermi Research Alliance, LLC under 
Contract No. De-AC02-07CH11359 with the United States Department of Energy. 

Appendix A: Vector Loop Calculation 

In this appendix, we present the explicit calculation of the vector loop contribution to gluon 
fusion. As the loop calculation involves a particular choice of gauge, we perform the calculation 
in both unitary (R^ — ?• oo) and Feynman {R^ = 1) gauges. This calculation generalizes the well- 
known Standard Model calculation for /i — t- 77 |45i [68l [83] to situations where a new "Higgs" 
field acquires a vev and leaves a non-Abelian gauge symmetry unbroken. In the SM, the Higgs 
field is responsible for spontaneously breaking SU{2)l x U{1)y, leaving the photon as the gauge 
field of the remaining Abelian U{l)em gauge symmetry. In contrast, in the renormalizable coloron 
model (ReCoM), <I> is responsible for spontaneously breaking SU{'i)i x 5[/(3)2, leaving the gluon 
as the gauge field of the remaining non-Abelian SU{3)c gauge symmetry. Then, gg — )• (pji is the 
non-Abelian mirror version of /i — >• 77. We intuit that M.{gg — )■ (pR) is enhanced by a color factor 
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over the mirror process A^(/i — )• 77), which is borne out from our calculation. 

1. Vector Loop Amplitude: Unitary Gauge 

We present the unitary gauge calculation of a colored vector contribution to gluon fusion. As 
mentioned above, we assume an extended color gauge symmetry that is partially broken by the 
vev of a new scalar field After Higgs portal-induced mixing of h and the matrix elements 
M.{gg —5- si) and M.{gg — )■ S2) are simply related to M.{gg — >■ h) and M.{gg — )■ </>_r) by 



M{gg si) = ce [M{gg /i)]L^=^^^ sg [M{gg (/'i?)]L,^=„ 
Migg S2) = se [Migg h)]\^^^^^^ +ce [M^gg L^^=„^ 



(Al) 



There are two diagrams which contribute to Ai{gg — )■ 0r) in the unitary gauge, shown in Fig. [7] 
of the main text. The triangle diagram in Fig. [7]A for the coloron insertion gives 

^Mf = -gl ( r'^f^'e,,e,. 



(27r)'^ (A;2-m2„)(Pa-"iG')(Pfe-"^G") 
where Pa = k + pi, pb = k — p2, and the three-point vector vertex is 

V>^-P = {k+paYg'P + {-Pa-PiYgP^ + {Pi - kYg^'^ 
^7/3" = + kYg'^P + _ p^Yg-yP + {-k - pifg'''' . 

The amplitude for the bubble loop in Fig. [Tfi is 



(A2) 



(A3) 



^■^^ - " [2) [-^) '^^'^^ J (2^ ipl-ml,Yk^-ml,) ' (^^^ 

where pa = P1+P2 — k, a symmetry factor of - has been included, and the four-point vector vertex 
is 

^bcyiJ-PlP _ _j^g2^bc ^ jace jdbe ^gfi-y gpl3 _ g^i/B gPf^ _|_ jade jcbe (^gUPg-yji _ gpl3g^p^ 

_^ jabejcde ^gPPgjfS _ gpjgpl3^^ (A5) 

= -iglr^^f^^{2g^'"<gp^-g^'^g'^-g^''g"<^) ■ 
After expanding both Eqs. ( |A2[ ) and ( |A4[ ) using Feynman parameters, performing the loop mo- 
mentum integration using Dimensional Regularization |84j . and some simpification, we arrive at 
the summed amplitude 

iM^v' = ^ ) (^) ^'''e^.e,, [pU - "^9^') Fyira^) , (A6) 
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where C{rQi) = 3 for the coloron, tg' = m'^_^ / (AuIq,) , and the loop function Fy is given in the 



main text in Eq. (55). We remark that this result also agrees with the analogous SM calculation 



for h — > 77 with the appropriate substitutions Og — )• a, C{r) — ?■ 1, — ?■ Vh, and a factor of 2 
w 

included for the charge conjugate process [l5l |68l |83] . 



2. Calculation: Feynman Gauge 



pi — •- ji. a 




Pi— 7,d 



(A) 




(B) 



(C) 




Gc 



(D) 




9 TrQT3T5T5Wtr 



(E) 




9 



9 THSTSISTSWO- 



(G) 



Pi — " ^. a 



9annroirsTnr*--'f 

P2 — - 7.ii 
(H) 




9 n5Tn5T575W5" 



(I) 




(J) 



FIG. 10. Feynman diagrams for vector loop contributions to gg — si in the Feynman gauge. 



As a cross-check of our result in Subsec. A 1 , we perform the same calculation in Feynman gauge, 



setting = 1. Again, we adopt the ReCoM with the Higgs portal, as discussed in Subsec. VIA 



In contrast to the calculation done in the unitary gauge, here we perform the calculation in the 
scalar mass basis, explicitly deriving the corresponding Goldstone couplings to si after taking into 
account the h-cpR mixing. This motivates an interesting discussion of the coupling of Goldstone 
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bosons to their partner Higgs field when the partner Higgs field is mixed with spectator scalars 
and provides a further check that Higgs mixing and NP contributions to si^2 production can be 



factored as in Eq. (Al). 



3. Goldstone couplings 



Beginning with the full scalar potential in Eq. (42) and the exact vacuum expectation values 
given in Eqs. (43) and (44), we examine the Goldstone couplings to the scalars h and (j). After 



expanding the potential, we have 
1 



(pR + y (3A<i. + fi:$)j + h{-XHVh)^ + ^ {-XhpV<j,(t>R + 2XHVhh) 



(A7) 



This can be written in matrix form as, 



h (t>R 
a/2 V2 



h (jrR 
^/2 ^/2 



( 2\hi 



p2 



M 



^0 



(A8) 



VhV2 



Note if we set X^p = 0, then Vh and v^j, become the unperturbed vevs and the Goldstone couplings 



become mj^/vh and m^^/v^ for the original m^^, mj^, v^, and Vh- From Eq. (A8) and Eq. (47) we see 



that the Goldstone-Goldstone-scalar interaction matrix M is the same as the scalar mass matrix. 
Thus when we diagonalize the mass matrix, we will simultaneously diagonalize the Goldstone 



interaction matrix in Eq. (A8). Explicitly, we write the scalar mass and Goldstone interaction 
terms as 




h 4>R 
V2 V2 



( 



M 



\ 



After diagonalization we obtain, for the interaction term. 



ri2 



V v^V2 J 
\ 



(A9) 



(AlO) 



V v^V2 J 

where Md = U~^MU is the diagonalized mass matrix with eigenvalues given in Eq. (50), and 



U is the unitary rotation matrix with its mixing angle defined in Eq. (49). We see that the 



32 



Goldstone-Goldsone-si interaction is 

2 

C D -se^siChGU"^ , (All) 

and so the Goldstone-Goldstone coupling to the scalar mass eigenstate is proportional to the scalar 
mass squared and a mixing angle. This justifies our ability to factor out the Higgs mixing angle 
when considering si^2 production. 



4. Continuation of the Feynman Gauge Calculation 



Returning to the calculation, there are ten diagrams in the Feynman gauge that must be cal- 



culated, as shown in Fig. 10 Using Pa = k + pi and pi^ = k — p2, Pc = Pi + P2 — k, Eq. (A3), 



and Eq. ( A5 ) , we obtain the amplitudes 



J ^ J {2tt)'^ [pI - ml,) (P - ml;,) {pi - m^,) 
^2 ^ d'^k 1 



--29./ / (^-».— j .1..,, y _ _ (A21) 

After converting to Feynman parameters, calculating the loop integrals, summing all of the ampli- 
tudes, and some simplification, the divergences cancel and we obtain for the production of si, 

iM^y^ = -^se (^) 6^'^e,,e2, [p^ " '^9^'') Mrc) , (A22) 
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in agreement with Eq. (A6) and Eq. (Al). 
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